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ABSTRACT 


A relaxation  process  is  applied  to  the  problem 
of  matching  relational  structures  involving  numeric 
quantities.  The  method  is  empirically  shown  to  be 
very  sensitive  to  the  exact  form  of  the  updating 
rule  used.  However,  with  the  proper  updating  rule, 
the  method  works  well,  and  is  remarkably  tolerant  of 
measurement  errors. 
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1.  Introduction 


i 


This  paper  addresses  the  following  problem:  Suppose  we 
have  a description  of  a certain  object  in  terms  of  its  parts, 
their  properties,  and  the  relations  between  them.  Now  suppose 
we  are  given  a description  of  some  part  of  the  real  world  in 
terms  of  these  same  properties  and  relations.  How  can  we 
determine  whether  an  instance  of  this  object  actually  occurs 
in  the  given  part  of  the  world?  A good  example  is  the  inter- 
pretation of  visual  images.  We  may  describe  a house,  for 
example,  in  terms  of  its  parts  (roof,  walls,  windows,  doors) , 
their  properties  (such  as  shape,  orientation,  color  and 
texture) , and  the  spatial  relations  between  them  (roof  above 
walls,  windows  surrounded  by  walls,  etc.).  Given  a visual 
scene,  we  may  segment  its  image  into  regions,  and  describe  the 
scene  in  terms  of  these  regions,  their  properties  and  their 
relations.  We  may  now  ask  the  question.  Is  there  a house  in 
this  scene?  That  is,  are  there  regions  in  the  scene  which  look 
like  roofs,  walls,  windows  etc.,  and  are  those  regions  in  the 
proper  spatial  relationships  to  form  a house? 

A previous  paper  (Kitchen  1978)  has  dealt  with  descriptions 
using  properties  and  relations  that  must  either  be  wholly  true 
or  wholly  false.  However,  these  are  not  fully  adequate  for 
describing  real-world  objects,  so  we  are  led  to  consider  prop- 
erties and  relations  which  can  take  on  numeric  values,  and 
along  with  this,  we  must  be  prepared  to  deal  with  errors  in 


such  values,  permitting  "approximate"  instances  of  a description 
of  an  object.  The  next  section  provides  a formalism  for  such 
"quantitative"  properties  and  relations,  and  indicates  how 
relaxation  methods  can  be  applied  to  matching  descriptions  in 
these  terms.  The  third  section  contains  an  empirical  study  of 
the  effectiveness  of  these  methods,  while  the  fourth  discusses 
the  results  of  this  study  and  ways  in  which  relaxation  methods 
can  be  further  extended. 


I 


! . Definitions  and  Description  of  Method 

The  following  provides  a formal  notation  for  descriptions 
of  an  object  (or  a scene)  in  terms  of  its  parts,  their  prop- 
erties, and  the  relations  between  them. 

Definition 

A quantitative  relational  structure  is  a triple  <X,#,<p> 
whose  components  are  as  follows: 

x is  a finite  set,  called  the  carrier  of  the  relational 
structure.  The  elements  of  X are  called  nodes , and  cor- 
respond to  the  atomic  parts  of  an  object  which  are  used 
in  its  description. 

$ is  a finite  set  of  quantities . Each  quantity  in  $ has 
associated  with  it  a pair  of  positive  integers,  called 
respectively  its  order  and  dimensionality,  indicated  by 
superscripts  in  parentheses.  Thus  Qv  ' ' fc<t>  indicates  that 
Q is  a quantity  in  with  order  m and  dimensionality  n. 

For  example,  a quantity  corresponding  to  the  ordinary 
Cartesian  coordinates  of  a point  would  have  order  1 and 
dimensionality  2,  since  it  assigns  to  single  points  pairs 
of  numbers.  On  the  other  hand,  a quantity  corresponding 
to  the  ordinary  Euclidean  distance  between  points  would  be 
of  order  2 and  dimensionality  1,  since  it  assigns  a single 
number  to  each  pair  of  points. 

^ is  a mapping  that  associates  with  each  a partial 


function  from  XM  to  Rn,  where  R is  the  set  of  real  numbers. 
Thus  determines  exactly  which  m-tuples  of  nodes  have 
values  associated  with  them  by  quantity  Q,  and  what  these 
values  are. 

Let  and  x^  , x2 , . . . , xmtX.  By  way  of  notation  we  write 

q (m, n)  (x^Xj, . . . ,xm)  to  indicate  the  value  of  at 

(x^  x2#...,xm),  provided  the  function  is  defined  at  this  point. 
Following  the  usual  convention,  when  no  confusion  can  arise,  a 
quantitative  relational  structure  will  be  denoted  merely  by 
the  name  of  its  carrier. 

Below,  the  symbol  M will  normally  be  used  to  denote  a quan- 
titative relational  structure  thought  of  as  a model,  that  is, 
a description  of  a certain  object  in  terms  of  its  parts  and 
their  properties.  Similarly,  the  symbol  W will  normally  be 
used  to  denote  a quantitative  relational  structure  which  des- 
cribes some  part  of  the  world. 

In  order  to  find  an  instance  of  an  object  in  some  part  of 
the  world,  it  is  necessary  to  match  up  the  parts  of  the  object 
with  the  corresponding  pieces  in  the  world.  However,  because 
of  errors  and  uncertainties  in  measurements,  we  must  be  pre- 
pared to  assign  different  degrees  of  confidence  to  different 
pairings  of  model  parts  with  world  parts.  It  therefore  becomes 
necessary  to  introduce  the  notions  of  fuzzy  truth  and  fuzzy 
logic.  In  ordinary  logic,  a given  proposition  can  be  assigned 
one  of  two  truth-values:  it  can  be  either  "true"  or  "false". 


The  usual  logical  connectives  such  as  "and"  and  "or"  can  be 
thought  of  as  truth-functions  which  take  as  arguments  either 
"true"  or  "false"*  and  return  either  "true"  or  "false"  as 
results,  according  to  the  truth-table  for  each.  In  the  most 
common  treatment  of  fuzzy  logic,  a given  proposition  can  have 
assigned  to  it  a real  number  in  the  range  zero  to  one,  reflect- 
ing the  degree  of  truth  in  that  proposition.  Similarly,  the 
truth-functions  are  extended  to  take  real-valued  arguments  and 
produce  real-valued  results.  Suitable  extensions  of  the  "and" 
function  are  xy  or  min  {x,y}  for  arguments  x and  y.  For  the 
"or"  function  we  can  use  max(x,^  or  1- ( 1-x) ( 1-y ) . For  a more 
detailed  treatment  the  reader  is  referred  to  Zadeh  (1965). 
Definition 

Let  M and  W be  quantitative  relational  structures.  An 
assignment  of  M to  W (the  order  being  material)  is  a function 
which  assigns  to  every  pair  (x,y)tM*W  a real  number  in  the 
range  zero  to  one.  Thus  an  assignment  indicates  the  degree  of 
confidence  we  have  in  each  pairing  of  model  node  with  world 
node. 

Definition 

Let  M and  W be  quantitative  relational  structures.  Let 
g(m,n)  be  a qUantity . Then  the  goodness  of  fit  function  for 
Q,  written  Yq*  is  a function  YQ:Rn><Rn~M0, 1 ] . If  utRn  is  a 
value  of  Q in  M,  and  vtRn  a value  of  Q in  W,  then  Yq(u,v) 


measures  how  closely  u and  v agree,  with  1.0  indicating 
perfect  agreement,  and  0.0  indicating  no  agreement  whatsoever. 
Clearly,  the  choice  of  is  determined  by  our  tolerance  for 
measurement  errors.  A suitable  form  for  Yq  is: 


where  and  v^  represent  the  ith  components  of  u and  v 
respectively,  and  si's2(***,sn  are  scalin9  parameters  associ- 
ated with  Q.  This  form  of  the  goodness  of  fit  function  will 
be  the  only  one  used  in  this  paper.  It  has  a peak  of  1.0  when 
u and  v are  exactly  equal,  and  falls  off  towards  zero  as  u and 
v move  apart. 

Definition 

Let  M and  W be  quantitative  relational  structures  with 
quantities  4>,  and  let  A be  an  assignment  of  M to  W.  Let 
x€M  and  yfcW.  Then  the  local  compatibility  of  x with  y,  with 
respect  to  A,  C(x,y,A)  is  given  by 


n 

xl'x2 

Q (m*n)  t<t> 
UtRn 

such  that 


x=x^  for  some  i 
Q (x^, • . . , xm ) =u 


U 

yx ,y2  * • • • ,ymtw 

, _n 
vfcR 

such  that 
y=y^  for  some  j 

Q<yi ym)=v 


(likimA(xk'yk))“YQ(u'v) 


where  0 and  U stand  for  the  operations  of  fuzzy  logical  "and" 
and  fuzzy  logical  "or"  respectively.  While  this  formula  may 
seem  rather  daunting  at  first  sight,  it  falls  naturally  into 
three  parts,  each  of  which  has  a rather  simple  justification. 

The  innermost  part 

<i«UA,v»k>lnV°'v> 

captures  the  following  notion:  Suppose  we  have  two  instances 
of  a quantity,  Q(x^ , . . . ,xm) »u  and  Qly^, . . . ,ym)=v.  Then  these 
two  instances  can  only  be  considered  "compatible"  if  their 
numeric  values  match  (hence  the  terms  y^(u,v))  and  the  assign- 
ment permits  all  the  pairings  (x^y^)  , (x2,y2)  • • • • * (xm,y m> 

(hence  the  term  ) . 

Moving  outwards,  the  fuzzy  logical  "or"  U,  embodies  the 
following:  Suppose  we  have  an  instance  of  a quantity  Q(x^, . . . ,xm) 
which  mentions  x.  Then  in  order  for  x to  be  "compatible"  with 
y,  there  must  be  in  the  world  a corresponding  instance 
Q(y^,...,ym)  which  mentions  y,  and  which  is  "compatible"  in 
the  sense  of  the  previous  paragraph. 

Lastly,  the  outermost  fuzzy  logical  "and"  0,  ensures  that 
this  "compatibility"  is  satisfied  for  every  quantity  instance 
in  the  model  which  mentions  x. 

Note  that  there  is  a different  local  compatibility  function 


for  each  choice  of  fuzzy  logical  truth  functions. 


I 


I 


3.  Relaxation 

With  the  above  ground-work,  we  can  define  a relaxation 
scheme  which  generates  a sequence  of  assignments 

between  two  quantitative  relational  structures  M and  W, 
using  the  following  recurrence: 

AQ(x,y)  = 1.0  for  all  xcM,  yCW 

A^(x,y)  = C(x,y,A^_j)  for  all  xtM,  yt-W,  i>0. 

That  is,  we  set  up  an  initial  assignment  Aq , which  permits 
all  pairings  of  model  nodes  with  world  nodes.  Then  we  use 
the  local  compatibility  of  each  pair  to  generate  a new  assign- 
ment. This  process  can  be  repeated  for  a fixed  number  of  steps, 
or  untij.  the  sequence  of  assignments  converges  in  some  suitable 
sense.  In  fact,  if  we  use  maximum  as  a fuzzy  logical  "or", 
and  minimum  as  a fuzzy  logical  "and"  in  the  local  compatibility 
function,  and  provided  we  assume  that  goodness  of  fit  functions 
are  always  greater  than  zero,  then  we  can  prove  these  two 
results  as  generalizations  of  those  in  Kitchen  (1978) : 

(1)  There  exists  j such  that  A.=A.  for  all  k>j. 

J K 

(2)  If  we  treat  the  M and  W as  discrete  relational 


structures  (that  is,  we  ignore  the  values  of  quant- 
ities and  regard  them  as  discrete  predicates)  then 
the  relation  R={ (x,y) | A^ (x,y) >0}  can  be  considered  as 


\ 


II 

V 

I: 

\ 


j 

,j 

I 


to  W. 


1 

. 

These  results  follow  easily  from  those  in  Kitchen  (1978) 
when  it  is  realized  that  local  compatibility  is  the  fuzzy 
analog  of  the  notion  weak  local  consistency  defined  therein. 

Under  conditions  other  than  those  stated  above,  no  such  re- 
sults can  be  simply  proved,  and  the  experimental  results 
described  below  suggest  that  no  such  results  in  fact  hold. 


4 . Implementation  and  Empirical  Results 

A computer  program  was  written  in  the  programming 
language  Pascal  to  implement  the  above  relaxation  scheme. 
However,  two  small  changes  were  made.  Firstly,  the  goodness 
of  fit  functions  were  modified  to  return  zero  if  any  term 
in  the  summation  were  greater  than  a certain  bound.  This  was 
necessary  to  prevent  arithmetic  overflow  when  comparing  widely 
disparate  quantity  values.  In  the  experiments  described  here, 
this  bound  was  set  at  100.  The  other  modification  permitted 
the  program  to  omit  from  further  processing  any  pairing  whose 
assignment  value  had  fallen  below  a certain  threshold.  This 
saved  computer  time  by  avoiding  re-computation  for  pairings 
with  so  little  support  that  they  could  be  disregarded.  The 
threshold  used  here  was  10  6 . In  several  trials,  the  use  of 
such  a threshold  had  a negligible  effect  on  the  course  of  the 
relaxation,  but  often  reduced  computation  time  by  a factor  of 
two  or  three. 

Several  experiments  were  run  to  investigate  the  behavior 
of  the  relaxation  under  various  conditions.  This  section 
describes  what  was  common  to  all  experiments. 

All  the  experiments  used  as  a world  a quantitative  rela- 
tional structure  describing  44  cities  of  the  eastern  United 
States  in  terms  of  their  populations  and  the  road  mileages 
between  them.  Thus  "population"  was  a quantity  of  order  1 


and  dimensionality  1.  which  was  defined  at  every  city, 
taking  values  between  0 and  7,605,000.  The  scaling  parameter 
s^  in  its  goodness  of  fit  function  was  10,000.  "Distance" 
was  a quantity  of  order  2 and  dimensionality  1,  defined  for 
every  pair  of  cities  with  a direct  road  route  between  them. 

It  took  values  between  37  and  387,  and  had  a scaling  parameter 
of  10.  The  data  were  drawn  from  the  Rand-McNally  1976  Road 
Atlas. 

Models  of  various  sizes  were  generated  by  randomly  select- 
ing compact,  connected  groups  of  cities,  and  then  including 
all  populations  of  these  cities  and  all  direct  distances  be- 
tween them.  A perturbation  of  size  d (say)  could  be  introduced 
into  a model  by  the  following  method.  For  every  quantity 
value  in  the  model,  generate  a random  number  uniformly  dis- 
tributed in  the  range  [-d,+d] , multiply  it  by  the  scaling  para- 
meter appropriate  to  that  quantity,  and  add  it  to  the  quantity's 
value.  Thus,  a perturbation  of  size  10  would  lead  to  an  average 
deviation  of  50  miles  in  distance  measurements,  which  is 
roughly  a 25%  relative  error.  This  feature  permits  investi- 
gation of  how  well  the  relaxation  copes  with  errors  in  meas- 
urement. 

The  first  experiment  tried  various  forms  of  the  local 
compatibility  function  under  increasing  perturbations  of  the 
models  used.  The  different  forms  were  as  follows: 


Form 

Is 

maximum 

for 

fuzzy 

logical 

"or" 

minimum 

for 

fuzzy 

logical 

"and" 

Form 

2: 

maximum 

for 

fuzzy 

logical 

"or" 

product 

for 

fuzzy 

logical 

"and" 

Form 

3: 

complemented 

-product  for 

fuzzy  logical 

product 

for 

fuzzy 

logical 

"and" 

The  complemented-product  has  the  formula 

l-(l-a)(l-b)  for  real  arguments  a and  b. 

Three  other  forms  (4,5,  and  6)  of  the  local  compatibility 
function  were  used.  These  corresponded  to  forms  1,  2 and  3, 
except  that  the  outermost  application  of  the  fuzzy  logical 
Mand"  in  the  local  compatibility  function  was  replaced  by  an 
arithmetic  mean.  Roughly  speaking,  these  forms  of  the  local 
compatibility  function  would  be  prepared  to  accept  evidence 
that  looks  good  on  average,  rather  than  insist  that  all  evi- 
dence be  found. 

Two  measurements  were  used  to  gauge  the  effectiveness  of 
the  relaxation.  The  first  is  called  the  "truth-value  difference". 
In  any  assignment,  there  will  be  some  pairings  which  are  cor- 
rect, and  many  pairings  which  are  incorrect.  Ideally,  all 
correct  pairings  should  receive  the  value  one,  all  incorrect 
pairings  zero.  This  ideal  situation  will  not  always  occur, 
but  a good  measure  of  how  well  an  assignment  separates  correct 
pairings  from  incorrect  pairings  can  be  obtained  by  taking 
the  average  truth-value  assigned  to  correct  pairings,  and 


! 


' 


r 


subtracting  from  it  the  average  truth-value  of  incorrect 
pairings.  However,  since  we  are  interested  in  the  separ- 
ation of  good  pairings  from  bad,  and  not  absolute  truth 
values,  the  values  used  for  computing  the  averages  were  not 
the  raw  truth-values  of  the  assignment;  instead  they 
were  linearly  scaled  by  division  by  the  largest  truth-value 

f 

in  the  assignment.  This  also  permitted  a more  equitable  comp- 

i 

arison  across  the  different  forms  of  local  compatibility, 
since  all  the  truth-values  were  standardized.  Note  that  this 
normalization  occurs  only  when  actually  computing  the  truth- 
value  difference;  the  assignment  used  for  the  relaxation  was 
in  no  way  altered.  In  the  best  case  the  truth-value  difference 
will  have  value  1.0,  at  very  worst  it  could  be  -1.0,  however 
in  the  normal  worst  case  it  would  take  a value  very  close  to 
zero,  indicating  that  the  assignment  makes  no  significant 
distinction  between  correct  and  incorrect  pairings. 

The  other  measurement  used  was  called  "good  separation". 

An  assignment  is  said  to  have  good  separation  when  the  correct 

I 

pairings  have  higher  truth  value  than  any  incorrect  pairing. 

Since  good  separation  is  either  present  or  absent,  it  cannot 
distinguish  between  those  assignments  which  only  just  fail  to 
have  good  separation  and  those  which  fail  miserably.  However 
it  captures  the  important  desideratum  that  the  correct  pairing 
of  model  with  world  can  be  found  merely  by  selecting  the  best 
pairs  in  the  assignment. 


For  each  form  of  the  local  compatibility  function,  four 
randomly  generated  models  of  size  5 were  used,  and  these 
models  were  perturbed  by  amounts  0 through  5,  and  10.  The 
relaxation  process  was  run  for  five  iterations,  and  measure- 
ments made  on  the  final  assignment.  Figure  1 shows  the  average 
truth-value  difference  over  the  four  models,  while  Figure  2 
shows  the  number  of  models  (out  of  the  four)  which  exhibited 
good  separation.  The  same  series  of  random  models  and  per- 
turbations was  used  for  each  form  of  the  local  compatibility 
function,  so  that  the  results  are  comparable  across  the  table 
(in  each  row,  all  six  forms  were  applied  to  the  very  same 
matching  problems).  Figures  3,  4 and  5 apply  only  to  form  1 
of  the  local  compatibility  function,  the  best  behaved  of  the 
six  tried. 

Figure  3 is  derived  from  the  same  experiment  as  Figures 
1 and  2.  It  shows  the  average  truth-value  difference  at  the 
end  of  each  iteration.  (The  rightmost  column  of  Figure  3 is 
identical  with  the  leftmost  column  of  Figure  1.)  Each  row  of 
Figure  3 gives  some  idea  of  the  convergence  behavior  of  the 
relaxation  process. 

Figures  4 and  5 demonstrate  the  effect  of  changing  the 
model  size.  For  each  of  three  sizes  5,  7,  and  10,  four  models 
were  generated  of  that  size.  Each  model  was  perturbed  by 
various  amounts.  Five  iterations  of  the  relaxation  process 


1 


were  applied,  and  measurements  were  taken  on  the  final 
assignment.  Figure  4 shows  the  truth-value  differences 
(averaged  over  the  four  models) . Figure  5 shows  the  number 
of  models  (out  of  four)  which  exhibited  good  separation. 
(The  leftmost  column  of  Figure  4 is  identical  with  that  of 
Figure  1;  similarly  for  Figures  5 and  2.) 


5.  Discussion  of  Results  and  Extensions 


One  of  the  most  striking  things  about  the  results  shown 
in  Figures  1 and  2 is  that  the  ability  of  the  relaxation 
process  to  tolerate  noise  depends  very  much  on  the  particular 
form  used  for  the  local  compatibility  function.  Those  forms 
which  use  multiplication  as  a fuzzy  logical  "and"  seem  to  fare 
the  worst.  While  remarkable  at  first  sight,  this  behavior 
has  a rather  simple  explanation,  which  is  best  illustrated 
by  an  example:  Suppose  we  combine  three  truth-values:  0.1, 

0.1  and  a,  where  a<0.1.  The  minimum  of  these  three  numbers 
will  reflect  directly  any  changes  in  a;  while  with  the  product, 
any  changes  in  a will  be  attenuated  by  a factor  of  0.1x0.1=0.01. 
Thus,  the  use  of  multiplication  will  not  produce  such  a good 
separation  of  truth-values  in  an  assignment,  even  when  a linear 
rescaling  is  done.  In  practice  it  is  observed  that  while  the 
truth  values  of  incorrect  pairings  are  reduced  quickly  to  zero, 
the  truth  values  of  pairings  without  perfect  support  (but 
which  are  nonetheless  almost  correct)  are  also  reduced  very 
quickly  to  zero.  However,  when  using  minimum  as  a fuzzy  logi- 
cal "and",  a pairing  without  perfect  support  is  reduced  to  a 
non-zero  level  which  reflects  the  strength  of  its  weakest  sup- 
port, and  stabilizes  at  that  level.  Thus  multiplication  works 
badly  because  it  too  readily  throws  out  pairings  which  should 


have  been  retained. 


On  the  other  hand,  the  use  of  an  average  at  the  outer- 
most level  of  the  local  compatibility  function  does  not  work 
as  well  as  the  use  of  a minimum  there,  because  the  relaxation 
then  holds  onto  incorrect  pairings  which  should  have  been 
discarded.  It  does  this  because  some  incorrect  pairings  may 
have  generally  good  support,  and  only  fail  in  one  or  two 
respects.  To  an  average  such  pairs  would  look  good,  while  a 
minimum  would  reject  them.  Note  however  that  these  remarks 
apply  only  to  the  experiments  described  above,  where  "noise" 
was  introduced  into  a model  by  perturbing  the  values  of  quan- 
tities without  otherwise  altering  the  structure.  It  is  also 
possible  to  introduce  "noise"  by  deleting  nodes  and  quantity 
instances  from  the  world.  In  visual  terms,  this  corresponds 
to  identifying  an  object,  parts  of  which  may  be  obscured  or 
otherwise  invisible.  For  this  type  of  problem,  the  averaging 
of  evidence  may  very  well  be  the  better  method,  as  is  suggested 
by  the  work  of  Ranade  and  Rosenfeld  (1978). 

Barrow  and  Tenenbaum  (1976)  in  their  paper  about  MSYS  remark 
that  the  particular  choice  of  fuzzy  logic  functions  seems  not 
to  have  much  effect  on  the  relaxation.  The  different  behavior 
observed  here  can  probably  be  attributed  to  the  fact  that 
these  experiments  deal  with  larger  matching  problems,  and  with 
apparently  larger  perturbations  than  theirs. 

In  Figure  3 we  see  the  improvement  in  the  separation  in  the 
assignment.  Not  only  does  perturbation  of  the  model  degrade 


the  result  at  a given  iteration,  but  it  also  increases  the 
number  of  iterations  required  to  reach  convergence. 

Figures  4 and  5 show  that  the  size  of  the  model  has 
little  effect  on  how  good  an  answer  the  relaxation  produces. 

It  should  be  remarked  in  passing  that  the  computation  time 
varies  in  a roughly  linear  fashion  with  the  size  of  the  model. 
This  is  because  the  computation  associated  with  each  model  node 
on  a given  iteration  is  fixed,  depending  only  on  the  quantity 
instances  which  explicitly  mention  that  node;  while  the 
number  of  iterations  required  for  convergence  seems  unaffected 
by  model  size. 

Relational  structures  can  be  extended  in  many  ways  to  make 
them  more  useful  tools.  Some  of  these  have  already  been 
mentioned  (Kitchen  1978).  However,  one  rather  simple  extension 
to  quantitative  relational  structures  would  seem  to  be  quite 
powerful.  First  of  all,  instead  of  having  a general  goodness 
of  fit  function  for  a quantity,  we  should  permit  a different 
goodness  of  fit  function  for  each  instance  of  a property  in 
a model.  By  this  means  we  could  enforce  a tighter  match  to 
more  critical  parts  of  a model,  while  easing  the  constraints 
on  less  critical  parts.  We  could  then  make  queries  like  "Find 
a city  with  exactly  68,000  people  about  100  miles  from  another 
city  with  a population  of  roughly  100,000."  Furthermore,  there 
is  no  need  for  goodness  of  fit  functions  to  be  symmetric  with 


L 


respect  to  deviations  above  and  below  the  ideal.  Using 
asymmetric  functions  we  could  pose  problems  such  as  "Find  a 
town  with  population  about  10,000  which  is  less  than  100  miles 
from  a city  of  at  least  500,000  people." 

Even  without  such  embellishments,  relaxation  seems  to  be 
a very  useful  method  for  matching  quantitative  relational 
structures  under  distortion.  With  the  use  of  a well-behaved 
local  compatibility  function,  the  method  is  remarkably  tolerant 
of  measurement  errors. 
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Figure  1.  Truth-value  difference  for  different  forms 
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